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I. 


OBJECTIVE 


The  objective  of  this  report  is  to  describe  data  collection 
and  computer  software  used  to  characterize  armored  targets  and  ground 
clutter  in  the  3-  to  5-nm  and  8-  to  14-nm  infrared  (IR)  spectrum.  The 
software  is  designed  to  evaluate  a conceptual  seeker  algorithm  for 
effectiveness  against  an  armored  threat  within  the  threat's  operational 
scenario. 


II.  INTRODUCTION 

The  US  Army  Missile  Research  and  Development  Command  (MIRADCOM) 
Advanced  Sensors  Directorate  is  conducting  a Ground  Target  Signatures 
Program  to  investigate  and  evaluate  passive  IR  targets  and  background 
(clutter)  signatures  and  identify  potential  target  discrimination  and 
acquisition  techniques.  To  these  ends,  several  computer  software  routines 
have  been  developed  which  progress  from  a point  of  producing  graphics 
which  provide  intuitive  insight  into  the  problem  to  qualitative  evalua- 
tion of  detection  and  acquisition  algorithms.  The  computer  software 
features  are  as  follows: 

a)  Three-dimensional  plot  of  target-clutter  energy  versus 
aspect  position  with  a subroutine  to  convert  the  data  to  a two-dimensional 
plot  of  energy  versus  scan  direction  position  — This  is  accomplished  by 

an  integration  process  in  the  vertical  dimension.  These  data  provide 
information  for  scan  direction  Weiner  Spectrum  analysis. 

b)  Two-dimensional  matrix  plot  of  energy  versus  target 
clutter  position  — This  program  provides  a matrix  of  data  values  versus 
position  to  be  used  in  a simple  hot  spot  detection  algorithm. 

c)  Target-clutter  noise  discrimination  model  used  to  vary  a 
threshold  on  the  matrix  array  provided  in  Routine  b)  — The  program  out- 
put is  the  number  of  data  points  exceeding  the  threshold  and  the  loca- 
tion of  those  points. 

d)  Histogram  of  target-clutter  scene.  This  routine  generates 
an  energy  probability  density  function  which  can  be  used  to  evaluate 
target  energy  versus  clutter  energy  and  the  setting  of  a seeker  threshold. 

e)  Variable  field-of-view  (FOV)  scan  — This  program  is  used 
to  subdivide  a large  array  of  data  into  smaller  increments  as  would  be 
seen  by  a matrix  array  of  detectors  or  a line  array  with  scan  direction 
data  storage  equivalent  to  the  vertical  scan  dimension. 

f)  Fly-in  simulation  model  — This  program  selects  consecutive 
subsets  of  the  original  data  matrix  to  simulate  the  reduced  FOV  encountered 
during  seeker  fly-in  and  calculates  data  statistics. 

g)  Graphics  target-clutter  frame-to-frame  comparison  model— 

— This  program  is  used  to  accumulate  statistics  over  a large  number  of 
scenes  as  in  Routine  d) . It  then  calculates  the  mean  and  standard 
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deviation  o£  the  accumulated  scenes  and  uses  Routine  c)  to  evaluate  a 
threshold  sat  to  the  scene  mean  value  plus  integer  multiples  of  the 
standard  deviation. 

h)  Two-dimensional  gradient  and  spatial  discriminator  — This 
program  generates  a two-dimensional  thermal  gradient  by  evaluating  the 
pixel  to  pixel  changes  in  both  horizontal  and  vertical  directions  through 
the  use  of  a 2-horizontal  by  2-vertical  pixel  window. 

i)  Two-dimensional  thermal  gradient  and  spatial  discriminator- 
This  routine  is  identical  to  Routine  h)  except  for  the  window  function 
which  has  been  expanded  to  3-horizontal  by  3-vertical  pixels. 

These  routines  are  used  in  conjunction  with  AGA  Thermovision  data 
to  evaluate  tactical  scenes.  The  following  sections  provide  the  AGA 
thermovision  raw  data  and  data  formating  into  a CDC  6600  compatible 
format,  a mathematical  and  physical  description  of  each  computer  program, 
and  a users  guide  for  each  program. 


III.  PROBLEM  DEFINITION 

A.  Introduction 

The  evaluation  of  seeker  performance  against  armored 
threats  requires  a thorough  understanding  of  the  seeker  algorithm  and 
the  nature  of  the  threat  and  its  expected  environment.  Once  this  basis 
is  established,  more  sophisticated  methods  of  target-background  dis- 
crimination can  be  applied. 

3.  System  Characterization 

The  first  step  which  must  be  taken  in  any  computer  evalua- 
tion is  characterization  of  the  system  to  be  analyzed.  To  do  this, 
either  the  target  must  be  characterized  by  some  model  or  actual  field 
measurements  must  be  obtained.  A model  is  usually  a tradeoff  between 
accurate  representation  of  the  target-background  and  the  mathematical 
complexity  which  the  analytic  system  can  handle.  Field  measurements 
are  real-world  situations  limited  only  by  the  similarity  of  the  data 
acquisition  system  to  possible  seekers  and  by  the  number  and  variety 
of  measurements  that  can  be  made.  Once  the  data  are  collected,  they 
can  be  analyzed  on  numerous  levels  from  intuitive  to  analytical. 


IV.  DATA  COLLECTION 

The  data  collection  system  is  a raster  type  scanner  manufactured 
in  Sweden  by  AGA  AKTIEBOLAG.  The  system  contains  two  major  subassemblies; 
the  camera  head  and  an  electronics  control/display  console.  The  camera 
head  consists  of  a silicon  lens  with  a 134-mm  focal  length  and  a maxi- 
mum aperture  opening  of  f/1.5.  The  optical  system  uses  a variable 


aperture  stop  to  control  the  FOV  and  an  image  plane  scanner  designed  to 
produce  a raster  scan  at  16  frames/sec  over  a 10°  by  10°  square  FOV. 

The  scan  rate  is  determined  by  the  280  vertical  line  2/1  interlace  raster 
with  140  unambiguous  data  points  per  horizontal  line  and  140  unambiguous 
vertical  line  resolution  capability.  The  single  detector  is  an  InSb 
photovoltaic  detector  operating  at  77°K  by  means  of  a nitrogen  dewar. 

The  detector  angular  subtense  for  both  horizontal  and  vertical  FOV  is 
1.3  mrad. 

The  detected  video  signal  is  supplied  to  the  video  amplifying  circuits 
which  amplify  and  filter  the  signal  to  drive  a small  cathode  ray  tube 
which  generates  a pictorial  result.  At  the  same  time,  the  video  signal 
is  supplied  to  an  analog-to-digital  converter.  The  analog-to-digital 
converter  presently  digitizes  approximately  one  frame  per  second  and 
stores  data  serially  in  a PCM  format  on  one  track  of  a 14-track  tape  on 
an  AMPEt-1300  tape  recorder.  Each  digitized  data  frame  consists  of  140 
vertical  lines  with  140  data  points  per  line  in  a 10-bit  word  for  each 
data  point.  This  14-track  output  tape,  with  one  track  of  serially 
packed  digitized  Thermovision  output  data  contains  140  X 140  X 10  bits 
of  information  for  each  frame  digitized.  This  information  is  then 
selected  on  a per  frame  basis  and  recorded  on  a digital  parallel  seven- 
track  tape  compatible  with  input  data  format  requirements  of  an  Army 
CDC-6600  digital  computer.  Further  and  more  detailed  descriptions  of 
recording  techniques  and  the  reduction  of  Thermovision  data  will  be 
addressed  in  the  ensuing  sections  of  this  report. 

The  AGA  Thermovision  system  is  relatively  small  and  may  be  mounted 
on  a helicopter,  elevated  tower,  or  installed  in  a fixture  at  ground 
level  to  view  a ground  target  or  background  scene.  In  support  of  many 
different  air  defense  and  ground  target  signature  applications,  the 
system  has  been  equipped  with  eight  different  bandpass  filters.  Each 
of  these  filters  represents  an  IR  bandpass  of  military  interest. 

A.  AGA  Thermovision  Field  Measurements 

In  general,  the  AGA  Thermovision  system  has  been  used 
in  IR  measurements  of  ground  targets  viewed  from  both  ground  level  posi- 
tions and  elevated  platforms.  This  system  of  data  collection  has  proven 
to  be  reliable  in  many  and  varied  applications.  For  example,  this  sys- 
tem has  been  used  by  the  Air  Force  to  record  in-flight  IR  signatures  of 
jet  aircraft  plumes.  The  Thermovision  system  was  mounted  in  a pod  on 
the  wing  of  a chase  plane. 

For  the  current  application,  the  Thermovisions  were  mounted  on  a heli- 
copter to  collect  ground  target  and  clutter  data  from  various  altitudes, 
ranges,  and  aspect  angles  near  vertical.  Due  to  the  interaction  between 
the  target  and  ground  clutter  and  a desire  to  gain  more  of  an  intuitive 
insight  into  the  problem  and  evaluate  seeker  concepts,  a large  quantity 
of  data  has  been  taken  for  ground  targets  and  clutter. 
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Data  of  these  types,  along  wi^h  computer  analysis,  can  evaluate  and 
establish  limits  on  seeker  systems  acquisition  techniques  which  operate 
predominately  on  energy  levels  and  spatial  frequency.  These  types  of 
data  and  computer  analyses  can  also  be  used  to  assess  the  feasibility 
of  automatic  target  cueing  technology  in  detecting  and  recognizing  tacti- 
cal targets  in  forward  looking  IR  (FLIR)  imaging  systems . 

B.  Data  Formatting 

The  purpose  of  this  section  is  to  outline  steps  required 
to  process  and  reduce  the  raw  Thermovision  data  to  a computer  compatible 
format.  The  block  diagram  of  Figure  1 will  be  used  to  describe  the 
required  process  from  the  point  of  data  inception  to  the  point  of  input- 
ting reduced  data  to  the  CDC-6600  for  analysis. 


STEP  1 STEP  2 STEP  3 


Figure  1.  AGA  Thermovision  data  recording  process. 


During  Step  1,  each  Thermovision  data  frame  is  composed  of  140  lines 
with  140  resolution  elements  per  line;  consequently,  if  each  data  point 
is  represented  with  a 10-bit  word  (16  X 140  X 140  X 10  bits  of  informa- 
tion), the  string  of  3.136  X 10^  bits/sec  would  require  recording  on 
magnetic  tape.  This  data  rate  is  well  beyond  the  tape  drive  operating 
capability;  therefore,  an  analog-to-digital  converter  was  developed  to 
digitize  one  frame  per  second,  approximately  196,000  bits  of  information 
per  frame.  This  requires  the  matching  of  lines  from  frame  to  frame  to 
reconstruct  one  frame  out  of  every  16  frames.  To  assure  correct  frame 
reconstruction,  the  last  two  data  points  at  the  end  of  each  frame  nave 
a special  coded  value.  Therefore,  at  the  end  of  Step  3,  a 14-track 
AMPEX-1300  tape  is  generated  with  digital  PCM  Thermovision  information 
on  only  one  track.  A physical  representation  of  the  one  track  is  shown 
in  Figure  2.  Each  line  is  composed  of  Data  Point  1 (Dl)  through  Data 
Point  140  (D140) ; each  data  point  is  represented  by  10  bits. 


FRAME  1 

LINE  NO.  1 LINE  140 


Dl 


D140,  Dl 


Figure  2.  Data  frame  format. 
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For  each  line.  Data  Points  139  and  140  have  the  same  special  bit 
pattern  to  denote  the  end  of  a line;  in  addition.  Data  Points  139  and 
140  of  Line  140  are  made  up  of  a spefcial  bit  pattern  to  identify  the  end 
of  a frame.  After  one  frame  is  recorded  on  the  tape,  a new  frame  is 
started;  this  process  repeats  itself  until  the  measurements  are  completed 
or  a tape  is  full. 

Steps  4,  5,  and  6 (Figure  3)  are  illustrated  to  show  how  the  data 
are  recorded  to  meet  CDC-6600  input  format  requirements.  During  Step  4, 
a decommutator  is  used  to  select  the  correct  coded  words  at  the  end  of 
each  line  and  at  the  end  of  each  frame  to  reconstruct  each  data  frame 
correctly.  Each  frame  of  data  is  input  to  a PDP-11  digital  computer 
(Step  5).  During  Step  6,  a PDP-11  digital  computer  takes  the  serial 
string  of  bits  and  generates  a seven-track  parallel  digital  packed  tape 
which  may  be  input  to  a CDC-6600  computer.  It  must  be  recalled  that  a 
10-bit  word  was  output  from  the  analog  to  digital  converter  to  represent 
each  10-bit  data  point.  This  10-bit  word  is  now  right-adjusted  in  each 
of  the  18-bit  words  output  from  the  PDP-11  computer.  Thus,  at  the  end 
of  Step  6,  a magnetic  digital  tape  composed  of  Thermovision  data  has 
been  generated  and  is  compatible  with  CDC-6600  software. 


STEP  * STEP  5 STEP  6 


Figure  3.  Data  conversion  to  CDC-6600  format. 


C.  Validation  of  Data  Format 

The  initial  task  in  the  Thermovision  data  analysis  is  to 
input  to  the  Eglin  BASES  program  the  digital  tape  generated  during  the 
previously  mentioned  six  steps  of  the  data  reduction  process.  The  10-bit 
input  data  are  now  contained  in  an  18-bit  word  and  must  be  selected, 
sorted,  and  shifted  into  a 60-bit  word.  After  processing  through  the 
BASES  program's  GETPIC  routine,  a new  seven-track  parallel  data  tape  is 
created,  which  is  made  up  of  60-bit  words,  with  five  12-bit  words  in 
each  60-bit  word.  Furthermore,  the  10  bits  of  meaningful  data  are 
packed  in  each  of  the  12  bits  of  information  as  in  the  example  of 
Figure  4. 


7 


i 


“i 


■i 


» 


XX  ► XX  ► XX  -« ► XX  -« ► XX  *4 

I 10  BIT  10  BIT  10  BIT  10  BIT  10  BIT  I 

Figure  4.  CDC  60-bit  word  data  format. 

After  using  the  GETPIC  routine  one  time  to  generate  the  new  60-bit 
word  data  tape,  the  BASES  program  may  be  run  any  number  of  times  with 
the  new  data  tape  without  using  the  GETPIC  routine  again.  Outputs  from 
the  BASES  program  applied  to  measured  data  are  tabulated  in  Table  1 for 
100  of  the  260  frames  of  data  contained  on  the  tape. 


V.  DATA  ANALYSIS 

The  simplest  level  at  which  data  can  be  analyzed  is  by  a two- 
dimensional  printout  of  data  counts  versus  the  two  spatial  coordinates 
(Figure  5) . This  gives  a quick  intuitive  feel  for  the  general  distri- 
bution of  energy  in  the  data  frame.  Maximum,  minimum,  and  rough  average 
energy  across  the  data  can  be  obtained.  Further  insight  can  be  gained 
by  plotting  energy  versus  position  using  a three-dimensional  plotting 
program  (Figure  6) . This  three-dimensional  plot  provides  a wealth  of 
qualitative  information.  From  the  three-dimensional  plot,  an  evalua- 
tion of  target  energy  as  compared  to  background  energy  can  be  obtained. 

A rough  idea  of  spatial  frequency  content  and  thermal  gradients  can 
also  be  determined.  The  three-dimensional  plot  also  quickly  reveals 
the  distribution  of  equal  energy  levels  (level  slicing). 

From  this  point,  there  are  basically  three  different  approaches  to 
target-background  discrimination:  spatial  filtering,  feature  selection- 
classification,  and  pattern  recognition  [ 1] . Spatial  filtering  can  be 
achieved  most  easily  by  scanning  with  an  array  of  detectors.  A linear 
array  of  detectors  whose  output  is  summed  while  being  scanned  at  right 
angles  across  a scene  will  produce  a one-dimensional  distribution  of 
energy  versus  scan  position  (Figure  6) . This  output  can  be  converted 
to  a Weiner  spectrum  using  the  Fourier  transform  [2],  The  Weiner  spec- 
trum can  then  be  used  to  characterize  target  and  background  signatures. 

The  importance  of  characterizing  IR  backgrounds  is  unquestionable, 
but  methods  for  doing  so  have  been  hotly  debated  [2-5].  Most  recent 
efforts  have  centered  around  modeling  IR  backgrounds  as  Gaussian  noise 
distributions.  Such  distributions  can  then  be  described  by  Weiner 
spectra,  autocorrelations,  or  line  scan  distributions . Thus,  the  main 
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thrust  of  data  acquisition  and  analysis  of  the  MIRADCOM  IR  ground  target 
signatures  program  has  been  to  confirm  or  refute  these  models. 

Feature  selection  and  classification  also  requires  knowledge  of  the 
statistical  energy  signatures  of  targets  and  backgrounds.  To  do  this, 
a histogram  of  the  number  of  data  points  versus  scene  intensity  must 
be  accumulated  (Figure  7) . From  the  histogram  (probability  density 
function) , target  and  background  signatures  can  be  fitted  to  statistical 
distributions,  usually  assumed  normal  or  Gaussian.  These  statistics 
will  then  define  error  probabilities  for  various  seeker-target-background 
combinations.  A qualitative  feel  for  the  performance  of  an  adaptive- 
threshold  seeker  algorithm  can  be  gained  by  thresholding  the  two- 
dimensional  matrix  plot  around  the  target  mean  at  integer  multiples  of 
the  target  standard  deviation  (Figure  8) . Further  characterization  of 
target  and  background  signatures  requires  that  statistics  be  accumulated 
on  a large  number  of  data  frames  to  increase  the  statistical  confidence 
in  the  signatures  obtained. 

The  first  step  in  performing  pattern  recognition  studies  is  to 
define  the  shapes  present  in  a particular  scene.  This  is  most  easily 
accomplished  by  spatial  differentiation,  which  essentially  enhances 
edges  of  objects  (this  is  also  useful  in  evaluating  the  performance  with 
edge  tracker  systems) . There  are  many  methods  for  edge  enhancement  but 
the  simplest  computationally  are  linear  matrix  approximations  to  the 
spatial  gradient.  These  methods  involve  using  either  a two  by  two  or 
three  by  three  matrix  as  a window  for  calculating  differences  in  abso- 
lute value  among  adjacent  data  points.  Matrix  windowing  also  allows 
simulation  of  various  configurations  of  shaped  detector  arrays.  By 
altering  the  size  of  the  matrix,  simulation  of  seeker  fly-in  can  also 
be  achieved. 


VI.  SOFTWARE  DESCRIPTIONS 

The  program  P2171  is  a versatile  package  which  provides  many 
alternate  methods  to  display  the  information  in  a particular  scene. 

The  program  provides  options  to  enable  the  user  to  read  packed  or 
unpacked  data.  Various  operational  modes  can  be  selected  due  to  the 
compartmentalized  structure  of  the  program.  Energy  levels  can  be  output 

-2  -1 

either  uncalibrated  as  counts,  or  calibrated  in  W cm  str  by  using 
previously  determined  calibration  constants  and  various  parameters  read 
from  the  input  data  tape.  Coordinates  can  be  output  in  three  modes: 
uncalibrated  (row  and  column) , calibrated  linear  (ft) , or  calibrated 
angular  (deg) . Provision  is  made  in  this  program  for  editing  the  data 
to  decrease  the  effects  of  noise  spikes  or  delete  obviously  bad  data. 

The  program  can  output  a summary  (Figure  9)  of  each  frame  of  data  which 
includes  the  following  data:  time  before  and  after  frame,  time  differ- 
ence, frame  number,  instrument  settings  (filter,  aperture,  gain,  etc.) 
maximum  and  average  counts,  location  of  maximum  count,  and  other 
internally  related  data. 
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Figure  9.  Data  summary 


Data  can  also  be  output  as  a two-dimensional  matrix  of  energy 
values.  Energy  and  position  units  can  be  selected  as  described  previously. 
Highlights  can  be  added  such  as  marking  maximum  count  or  blanking  out 
cells  with  counts  below  a selectable  threshold.  The  matrix  plot  also 
provides  the  sum  of  all  counts  in  the  frame.  A portion  of  the  program 
provides  the  two-dimensional  matrix  output  for  the  first  intuitive 
examination  of  the  data.  After  the  target  and  background  statistical 
distributions  have  been  determined,  the  two-dimensional  matrix  can  be 
thresholded  and  output  to  give  a qualitative  feeling  for  seeker  effective- 
ness and  probable  errors  (Figure  10). 


This  two-dimensional  matrix  of  energy  values  can  be  plotted  in 
three  dimensions,  that  is  a plot  of  energy  (as  the  vertical  coordinate) 
versus  spatial  coordinates  with  the  benefits  described  earlier.  The  sum 
of  the  columns  is  also  plotted  as  a two-dimensional  projection  along 
the  edge  of  the  three-dimensional  plot  to  provide  a simulation  of  scan 
direction  line  detectors. 


The  KPLOT  program  provides  plots  of  the  statistics  of  single  or 
accumulated  data  frames.  It  can  produce  bar  or  point  histogram  plots 
of  the  data  probability  densities,  as  well  as  the  mean,  variance,  and 
standard  deviation  of  a normal  distribution  which  best  fits  those  data. 

Program  GOWENS  is  a fly-in  simulation  model  and  edge-tracker  model. 
The  fly-in  simulation  is  implemented  by  iteratively  reducing  the  full 
frame  of  data  to  a smaller  subset  of  data  points  to  simulate  range 
closure,  simultaneously  finding  the  statistics  for  the  reduced  scene. 

The  same  program  can  simulate  an  edge  tracker  by  the  same  process  on  a 
line-by-line  or  column-by-column  basis.  Level  slicing  and  thresholding 
can  also  be  done  in  this  program  using  a subroutine  which  places  only 
those  data  points  with  counts  inside  a gate  into  a reduced  matrix. 
Coordinates  are  encoded  with  the  data  values  so  that  the  gated  data  may 
be  accurately  reconstructed. 

Program  Gradient  does  the  2x2  and  3X3  spatial  gradient  calcula- 
tions. The  program  utilizes  a moving  matrix  and  an  absolute  value  cal- 
culation to  approximate  the  temperature  gradient.  The  2X2  method 
calculates  from  the  matrix: 


a b 


c d 


The  gradient,  s,  at  point,  a,  is  then  given  by 
S = | a — d | -4-  | b — c | 

where  ||  is  absolute  value. 
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Figure  10.  Threshold  matrix  sample  plot. 


The  3X3  method  uses  the  matrix 


a b c 


c e f 


8 h i 


and  calculates  the  gradient,  s,  at  point  e as 

s ■ Is* I + Is J 


where 

Sx  = (c  + 2f  + i)  - (a  + 2d  + g) 

Sy  - (g  + 2h  + i)  - (a  + 2b  + c) 

and  approximate  the  partial  derivatives  in  the  x and  y directions, 
respectively,  and  can  thus  be  used  to  simulate  edge  trackers. 


VII.  SOFTWARE  USER  INFORMATION 


A.  Introduction 


This  portion  of  the  report  is  to  provide  a guide  for 
prospective  users  of  Program  P2171,  Program  GOWENS  and  Program  KPLOT. 

Setup  procedures  and  various  types  of  input  and  output  will  be  demonstrated. 
It  should  also  be  noted  that  Program  P2171  was  modified  from  the  original 
Eglin  version  to  simulate  IR-guided  missile  systems  against  certain 
armored  targets. 

B.  Program  P2171  Setup 

The  following  files  are  used: 


Files  Description 

1.  Input/Tape  5 Cards 

2.  Output/Tape  6 Time  history  listing 


3.  Tape  11 


4.  Tape  10 


Primary  data  source  PDP-15,  picture 
data  tapes  or  previous  compressed 
P2171 

P2385  or  P2204  trackpoint  range 
data  tape 


Restrictions 
Always 
Always 
A lways 


Optional 
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Files  Description  Restrictions 


5.  Tape  12  Compressed  output  tape  of  input  Optional 

Tape  11 

6.  FILM&L  Three-dimensional  picture  plots  Optional 

7.  Tape  9 Picture  matrix  listing  Optional 
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190  INPU1 


ALL  CAPPS  JCf  IDENTIFIED  B*  HOLLFPTTM  COOES  IN  COLUMNS  1-5, 

-TOVWVS — Tl^",5 — CONT  A I w-DA tm— FtrAnS-0gPg*>0IN6  BN  CAPO 

TYPE.  FIELDS  LFFT  BLANK  AND  CARPS  OMITTED  ASSUME  A DEFAULT 


VALUE. 


CAPT  TY°PS 

TITLE  CAPO  - INOICATES  TTTLF  PARC  TO  FOLLOW 
tOgP«ULT  PLANK  TTTfl 


_2 "LOT  CARO 


(S E F AULTPN0TPL0YT0U?°UT) 


_J PRINT  CARC CECINES  PRINT  OPTIONS 


(DEFAULT  NO  PRINT  OUTPUT) 


-V TAPE  CAPO CEFINE“?^TAPg  »A<AMFTfpS,  OUTPUT  TAP?  OPTION 


l,  -r  lIT'  ) I AA  t I ’j*  UUI  FIJI  IM-t 

CEFAULT  - NUMTfo  fA*<PLES/RECORO  = 5 
NUHRfo  wILLIFFCDNCS/SAMPtF  = 12.31)5 


Nfl  PACKEC  0UroU I TAPE 


TIME  CARO 


DEFINE6  START, STOP  TIME  TO  PROCESS 


(OEFAULT  ENTIRE  TAP F ) 


E LIMIT"  CAKU  - UEFINTS  LtNE  IU  PICK  AS  FIRST  LINE  OF  THERMOVISIUN  PICTUPE 

(DEFAULT  11) 


T CLASS  C 


ARC  - CEFINES  CLASSIFTCATION  OF  PLOT  OUTPUT 
(OFFAULT  UNCLASSIFIcO) 


TA  MSN  LA9eL  CARC  - ONLY  ON  CLASS  CARO  OPTION 


70  AUIMORITY  C AV  U - ONLY  ON  CLASS  CARO  OPTION 
7 C 03WNGR10F  CART  - ONI  Y PM  Cl  ASS  CARP  OPTION 


8 MATRIX  CARO-  CEFINES  MATRIX  norlCNS 
(OEFAULT  -NO  MATRIX* 


9 TPACX  CAPO  - CEFI NE  SOURCE  QF  RANGE.  RFL  BRNG  ANO/OR  TRACK  PT 


1 fl  '(At  TO  CARP  - DFFTUFS  VALTr  q&TA  SWITCH  SFTTTNC,. 


(DEFAULT-  IGNORE  VALID  DATA  SWITCH) 


11  CONV 6-APO RPFGIFIE S--GONVER  f TON  ■ FAC  TOR  COR  A GIVEN  FILTER, 


GAIN,  ANC  APE  R ATtl°E  , 
cefault-  conversion  FACTOR 


12  DELETE 


CARO-  SPECIFIES  nELETr  TIMES  ,/ 
(OFFAULT  NO  QF1  FTRS) 1_ 


13  CENTER  CAPO-  SPFCIFIES  LIN  ANT  WCRO  OF  00 °E SIGHT  OF  AGA. 


(OEFAULT  LINE=6P,  wrap  n(,5) 


14  CAL  ONL v 


SRECirlFS  CAL3°ATI0N  FRAMES  ONLY 


15  MIN  CARO 


RRECIFTES  Tmf  of  A X VALUE  BELOW  WHICH  A PICTURE  IS 
_CONSXO£RE.Q__Tfl  _?E  NOl-Of  SIGNIFICANT  LFVEL  AND 


WILL  NOT  «F  BL0T TF  0 
(OEFAULT. = 4*DACXCRPUND  LEVEL) 


15  EDIT  CARO 


CEFINFS  EDIT  OPTIONS  AND  BACKGROUND  LEVELS 
(OEFAULT-ffc  OFF  I N IT  ton  qf  THIS  CAPO) 


17  LABEL  CAPO  - CFFINES  LOCATION  ANO  HOLLERITH  LABEL  TO  BE  PUT  CN 
FLOT._ 


18  G3  C*°0- 


(OEF  (ULT  - BLANK  LABEL  FOR  UNDEFINED  LOCATIONS) 
-STSNALS  P °PG  RAM  TO-Rf&IN  EXF  GUT  ION 


IB  TAIL  CA°D 


- SPECIFTfs  A CDMRTanT  LOCATION  CF  TAILPIPE  ANC  TP»CK 
POINT  (CEFAULT,-  NO^CONSTANT  LOCATIONT 


TXTLC  CABO 


UTIF 


"THIS  CAPO— 

RE  PUT  ON 

iaragag«*!w  wAvg^g^Rvan" 

INFORHATION  TO 

POTTCARO 

FORMAT 

t-4  SARD 

-PLOT- 

12 

11-12  IPLOT 

C - NQ  PLOT  OyTPUT  (DEFAULT) 

| * 30  PLOT’S  *»LUS  HiX  SCAL€  SELECT 

Xejd  Lri5,7j^i^Fui_^!OHi_ — 

-»«U- 


wou 


•default  in 

■ HMF  M'JKafR  T r>  ST  OS  10  fiiOT  ■ 
(DEFAULT  110) 

N0R0  ' TO  "FGIN  3D  PLOT 

iO'Mtlll  1» 

3FD 


-12- 


17-18  1»T0*»1 

I I r 

19-20  IP30W2  W0P3  NUH0PO  TO  STOP  30  PLOT 
a — Ex^By-atMca  cpa 


12 


(OEFAUtT) 
-RAMS 


12 
~r t 


3 * EVEP7  TNIRO  FRAME 


23-24  IV  JEW 


29-39 

31-40 


ISUH 
NT  OP 


SPECIFIES  VIEWING  CIRECTION 

3 * VIEW  FROM  TOP  OF  »G» 

1FH  POOH  LFFT-gc*6* 
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1 = PLOTS  S'lM  IN  M/STR  ON  3-0  PLOT 

SCALE  INPUT  Ic  IFLOT  5T  1 

1N-P AST ER  TOUNTS  0 ,LE.  MTOP.LE.  1924 


12 

F10.0 


PRINT  CARO 


COLUMN  NAMF 
-1-5 KA-RO- 


11-12 


CESCRTPTION 
SPRINT-* 


FORMAT 
A5- 


IP9INT  0 = NO  PRINT  OUTPUT  (0EC AULT ) 12 

1 = PRINT  TIMES.  FILTER,  APERATURf,  GAIN,  CAL IRPATE, 

SOLID  OAIA,  LENGTHS,  ANO  A~~SPECIFIEO  cell 

•SEE  I PVL , IPV4  - COL  23-26)  FOP  EACH  PICTURE. 

7 = PRINT  OPTION  1 PLUS  POINT  CELLS  IN  SFEC1FLC  LINES 
ANO  WOODS  AT  A fofCIFIEO  RATS . «SEE  COL  13-22) 
IPPINT.NE.2.  <*OL  17-22  ARE  If 


IF  .....  _ ...  __ 

13 -14 TPRL1 UNe-N'twogP- TO  BEGIN  B8XMT 


IGNOPEC. 


15-15  IPPL2 


17-19 


TPRH1 


(OEEAULT  11) 

LINC  NUMBF®  TO  STOP  PRINT 
"toec9ult  nm 

WORT  NIJNRPP  TO  PFOln  P3INT 
ODEFAULT  1L 


13- 


19-90 IPRW2 WORT  NUMoro  TO  STOP  POINT 
•OEFAULT  100) 

-31—2-? NP°I-NT ROEG!  f 1 6S— P3TNT  JMC  . 


12 

12 

"17 

-42 


1 = FVF4V  FRAwF  (OEFAULT) 

2 * E V c ft  otwco  P041.C 

~T~*  EVERT  thi-»<j-pi***£ 


73-24 


IP  VL 


LINC  NUMBC®  riE  PARTICULAR  CELL  TO  BE  PRINTEO  • 
I0ECAIIL  T CC1 

-W WT-NUxncc-oc  PARTICULAR  CELL  TO  8E  PRI».T60.- 
(OEFAULT  5 P) 


u 


~T  APE  ~CARO 


COLUMN 

ii:«_ 

J S--/A 


..NAME DESCRIPTION. 


KAPO  TAPE 
IT  APE 


.FORMAT 


Ja=EN0  PACKFO_tipc  OUTPUT  «nECAULT) ,T/P  ; , F"PlS  T*pr  Mr 

? * I/F  IS  PACSFp  TA°E 

NUMOFP  EAMPtFS/PFCTRO  ON  INPUT  TAPE ** 

ITSENO  /l^MAKE TFXT£ONAL  TA»E  TO  BE  SENT  FROM  TS*  <SrE  1%. 

Pwt  r rn  a ■ ▼•nr  n ir ) . . ..  ■ - - 


NSR 


o»ie 


CFLTAT t IMPL9ETWFFN/S»M0LES  ON  input  TAPE  IN 
MILLISECONDS.  (PFFAIJLT  12.305  MS) 


1 


20 


COLUMN 
1 -4 


DESCRIPTION 

LIMIT" 


( OEcAUL  T 11> 


COL'IMN  NA  ME 


11-12  ICLASS 


CESCRT  PTI  ON . ■ ' 


UNCLASSIcIcO  I DEFAULT) 


2 * SECRET 

3 = TOP  SPCRpt 


2-5=NUM°FP  OP  ED  PC  I AL  DOWNf.RAOF  CAPO  IN°UT  ♦ 
ll-l«3=NUMgPR  OP  PUR c SPECIAL  OOHNORAOE  CAROS 


FORMAT 

A5 


FORMAT 


M J 
12 


authority  caro 

Tm'IST  FOLLOW  MSN~l7peL  rA°0  IF  1°  = 1 THRU  5) 


-COLUMN NAME  CFSCRIotIO* 

1-59  IAIJTHR  CLASSTFIHG  AUTHOPITY  FOR  OOWNGRACE  STAMP 

_7 0^73 J1CLAS CECLASSIFTCATTON  OAJE  

80  M C = IF  PLANK  IOCLAS  ISTJSEO 

L = UPON  ' C TIFTC  A T TON  BY  THE  ORIGINATOR 
k --POT  AUTOMATICALLY -rrCLASSlFIEQ 


-pflPMAY- 
. AS9 

AlO_ 

Al 


OOWNGRAOE  CARO 

(MUST_FOLLOW  AUTHORITY  CARC  IF*  IR»2,  MSN  LABEL  CARO  IP  IR>16)  " 

■COLUMN NAME  -OX&CglJ  W ECRU AT 

1-80  HOLLFR I T H OCWNG°AOE  STAMP  TO  BE  USEO  IF  IR  > 10  8A10 

OR  KOntFICATTON  TO  OOWNGRAOF  STAMP  STARTING  WITH 

fKCtrASlF  1 0 A T I CNtrTNF  »-tr  C A RUS-UNtE SS-XCtr^TC— 9E 

MOD  THEN  6 CAROS  AL  LONE  0 IF  IR  >10,9  CAROS  ALLOWED 


COLUMN 

NAME 

OFSCRIPTIOM 

FORMAT 

1 -4 
11-12 

IMTX 

M ATR 

=0,  INDICATES  NO  MAtRIX  OUTPUT  (DEFAULT) 

»N,  IUOICATFg  FFLFCT  EVFRY  NTM  FRAME  THAT  TS 

Atf 

12 

13-14 

Ic  EET 

CONSIOf RFO  FOP  OUTPUT  (SEE  COL  15-16) 

=0,  «ATPIX  KTLLonr  OUTPUT  IN^OEG REE  CELL  UNITS 

12 

15-16 

RBCELL 

CONS I Jc P FOR  mat  d i x OUTPUT  IF'RELATIVE  BEARING 
HAS  CHANGFO  BV  Mnog  THAN  R9CELL  (OEG) . NOTE  1 

12 

17-18 

19  -?n 

IFINO  1 
TPRMTX 

r ipROGPAM,,SECEOTSr  OfcN,'Ma  f RIX  LE NG  TM  lpAR A “FTFRS 

1 i MATRTX  lTRPOFTTTON  T(1  TAPF9 

12 

TP 

31-40 

HLEPT 

IF  NE  0 , T HrN  IS  left  LI«II  OF  MATRIX  EFOM  THE 
TAIL  PIPE  LOCATION  (OCFAULT  = 2.0) 

no. z 

5}  -60 

NF  RMT  X 

FROM  TAIL  °rPE  LOCATTOV  (DEFAULT  = 8.8) 

IF.NE.G  THEM  TS.HAXTMIJm  NO. OF  FRAMES  TO 

F10.2 

Uit  r u 1“  to  pct,  i u*  hi 

NOTE!  POOOPJM  WT|  l S’Act  AT  CFLL  C ORRE SPOMOI NG 

TO  OLEFT  A N3  CONS  TRUST  TWF  mATRTX  TC  THE  CELL 

COROESPONOING  TO  ORIGHr  OR  NFRNTX  PRAmEA 

TRACK  CARO 

FORMAT 

COLUMN 
1 -4 

NAMF 

<»’0 

DESCRIPTION 
‘•TRAC1* 

A4 

11-1? 

IT  RK 

=-?,  OO  NOT  pf AO  TAFEin.  USE  MAX  A'CA  VALUE  T6 
FIND  TAIL  rT°F  LOCATION 

li 

EOF  TPACK  OOTNT  AND  T AI C"  PIPE  ( $F  F ' TA  I L CAPO) 

= 0,  oo  NOT  RrA n TA  = C1P.  OO  NOT  OUTPUT  IP  MATRIX 

= 1 , — RE  An— TAPEl"^  TAIL  PT'®? — tOCATIUN— TO — ng— SFLECTCO 

FROM  TV  VTOgO  LOCATION.  WHEN  NOT  AVAILABLE, 

USE  Mtx  AG A LOCATION.  (T=ACK  »T  from  VCO  SIGNAL! 

= 2,  RF  AO  T AP  - IP  , IIS  F MAX  AGA  TO  FINO  TAIL  PTPE, 
TV  V I OF  0 whfn  LOCATION  F = OH  m»y  AGA  CANNOT 

US' 

°E  found. 

13-14 

I TP  T 

1 i TRACK  POINT  HARK  to  gE  TAKEN  FPQM  T V VIDEO 

12 

22 


COLUMN 

NAME 

OESCRI°TI ON 

FORMAT 

r-R — 

11-12 

T .1  -i  (. 

ti 

“CflNV“  " 

FILTE®  nijmb-o 

APFOATUFF  niihRFP 

12 

12 

16  -lb 
31-40 

RJjflj 

GAIN  N1IMOFR 

CONVE ° S I ON  FACTO®  from  COUNTS  TO  N1(CM*»2»SR» 

cno  r r 11  T r. rrina  t»i  itthu 

I? 

*10  • 2 

NOTEt  roNVFusiON  FACTOR  EOUALS 

(AVERAGE  INTERNAL  L A3  COUNT )«< SLOPE )/( A VG 

OELETE  CARO 


COLUMN  CESCRI°TTON 

KARO  "OELF" 


21-29  T3ELE  PHMMFSSSS  = £>10  OF  C EL  F T E INTERVAL 


FORMAT 
A A 


I .IIIUMIIj 


CENTER  CAOO 


T N C c N T W0»0  NUMRFC  oe  CFNTfr  OF  AGA  ( OFF AUL T- A 5 I 12 

RFRLI*  RESOLUTION  OF  ONF  LINE  OF  AGA  ( OFF  AtJLT  = 1 0/ 69  OEG>  fiq.2 


CAL  ONLY  FARO 


NAMF  CESCRTPTIPN 

KARO  "CAL" 


MIN  CARO 


FORMAT 
4 


COLUMN  NAMF 

1-4 *4*0-- 

11-12  IE07 

— tt-ri  iron0 


OESORTPTION  FORMAT 

-“EOIT“ — Air— 

-0,  00  MOT  EOTT  NOISE  SPIKES  IN  PICTURE  12 

=1,  FOI T NOISE  SPIKES  <OEFAULt) 

'IE  GT  l.C,  THEN  SPECIfIES  COUNT  VALUE- ABOVE- WHICH- T1D  ,-2 
A CELL  WILL  BE  '*o»'SIOERBO  POSSIBLE  NOISE  SPIKE. 

- MAx<<4.f  Al/G.-Of  aAC*SPQUNO-P£  AKS)  , Jt)  ) 

HEN  SPECIFIES  VALUE  THAT  THF  F10.2 

ADJACENT  LINES  MUST  9 E 9FL0W 
"E  CONSIDEREO  A-  NOISE -SPIKE 
OUNO) 

COUNT. VALUE  IQ  BE 
EVEL 

0 OETERMINE  AND 


NOTE!  LET  AVC=  AVERAGE  OF  PF4KS  IN  FJC TURE  BELCH 
20"  COUNTS 

THE  OP-T-IHAL — pACKGR0UN3=  A V6  ♦ PCT*  AVC 

(OEFAULT  FOP  OCT  =.75) 


COLUMN  NAME 


11-12 
44  4 


CESCPI °TION 

-“CA9E“ : 

LINE  NO  FOO  PLACEMENT  Of  THIS  L«°EL  ON  30  PLOTS 

>OLLE*TTH  LABEL-I.NCQP“ATION  TO  PUT  AT  I LAB 

LOCATION  NOTc I up  TO  25  LABEL  CAPOS  MAY  BF  INPUT. 
AFTER  A LABEL  CAPO  HAS  BEEN  INPUT, THEN  IT 


FORMAT 


AZIMUTH  ANO  ELEVATION  Oc  TAIL  PI pp  AND  TRACK  POINT 
HERE  KNOWN  FO°  A PARTICULAR  PICTURE 


TAPFlt  (TNPUT) 


THF  PROG® AM  WILL  ACOPPT  QMC  OF  TWO  Types  OF  TAP?  AS  INPUT  ON 
DtIT-tT  rrSE^“TAP£ — CAROn 

TYac  1.  ( WHCN  LTAPC=,3* 

TThTS  I S A” T A F Fi °R  JO  U r FTC  a y~p6Trr,T4M  6<!20C  ON  F»  0 P 1 <5  Cp^PUTEPT 


NOPMALL  THIS  TAPP  WILL 
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C . Program  GOWENS  Setup 


This  program  was  developed  in-house  as  a simulation  tool 
to  evaluate  each  frame  of  IR  data  efficiently  without  utilizing  large 
amounts  of  computer  resources  (memory,  time,  etc.)  which  are  required 


by  Program  P2171. 
as  Program  P2171; 
Program  KPLOT. 

Program  GOWENS  will  read  the  same  input  tape  (Tape  11) 
sort  and  process  IR  data  are  to  be  plotted  with 

The  following  files  are  used: 

Files 

Description 

Restriction 

1. 

Input/Tape  5 

Cards 

Always 

2. 

Input/Tape  1 

Primary  data  source,  same  as  Tape  11, 
P2171 

Always 

3. 

Output/Tape  6 

Listing  of  variables  (user  may  select 
any  variable) 

Optional 

4. 

Tape  2/Tape  3 

Sort  merge 

Always 

5. 

Tape  7 

Processed  IR  data  (cataloged) 

Always 

Program  GOWENS  Input/Output 

1. 

KT  Card  . . . 
be  evaluated. 

Frame  number  card,  1615,  to  identify  those  frames  to 

If  frame  = 0 or  -1,  program  will  stop. 

2. 

IFIL  Card  . . 
time,  filter, 

. Time  ID  card,  5I3,2F10.3,  4A10,  select 
gain,  and  aperture  setting.  The  view  is 

frame  by 
also  included. 

3.  Tape  1 . . . Primary  IR  data  source  (Eglin  AFB)  . 

4.  Tape  2 . . . File  to  be  sorted,  input  to  Tape  3. 

5.  Tape  3 . . . File  containing  sorted  IR  data  (maximal  to  minimal 

intensity) . 

6.  Tape  7 . . . Output  of  processed  IR  data,  contains  frame  times,  the 

view  if  any,  largest  200  data  values,  and  the  location  of  each 
value  in  the  frame. 


D.  Program  KPLOT 

Program  KPLOT  is  a mathematical  model  developed  for  graph- 
ical evaluation  of  the  output  produced  by  Program  GOWENS . KPLOT  generates 
bar  and  point  plots  by  interacting  with  the  Tektronix  software  plotting 
routines.  In  bar  plot  graphing,  the  number  of  IR  data  points  is  plotted 


versus  the  intensity  range  which  can  be  expressed  in  measurements  of 
counts  or  watts  per  steradian.  KPLOT  will  also  illustrate  a frame  of 
data  (100  X 100  matrix)  and  show  the  location  and  intensity  of  all  hot 
spots  in  the  frame. 

The  following  files  are  used: 


Files 

1.  Input/Tape  5 

2.  Input/Tape  7 

3.  Input/AGII 

4.  Output/Tape  6 


Description 

Interaction  with  terminal 
Processed  IR  data 
Tektronix  Software 
Information  display  (Plots) 


Restrictions 

Always 

Always 

Always 

Always 


Program  KPLOT  Input/Output 


NOTE  KPLOT  is  an  interactive  program  that  is  run  on  a terminal. 

1.  Tape  7 . . . Primary  input  for  KPLOT  (same  as  Tape  7 in  GOWENS). 

2.  AGII  . . . Tektronix  Software. 

3.  Tape  5/Tape  6 . . . interactive  I/O  with  terminal. 

4.  Output  . . . barplots,  display  of  frame  and  location  of  heat  sources. 
(Figures  1 and  2) . 


VIII.  COMPUTER  PROGRAM  LISTINGS 

The  listing  for  BASES  Program  P2171  is  presented  in  Appendix  A. 
Program  GOWENS  follows  in  Appendix  B.  Program  KPLOT  listings  are  con- 
tained in  Appendix  C. 


IX.  SUMMARY 

The  nine  different  computer  models  outlined  and  described  in 
the  introduction  all  utilize  different  discriminants  to  investigate  the 
spatial  radiance  of  target  and  clutter.  Of  the  nine  different  computer 
models,  only  three  are  contained  in  the  Computer  Program  Listing 
(Appendices  A,  B,  and  C)  in  order  to  minimize  the  size  of  this  report. 

In  the  case  of  energy  computer  model  described  in  the  report  the  over- 
lapping areas  between  target  and  clutter  represent  a loss  function  where 
targets  would  be  classified  as  clutter  or  vice  versa.  In  the  search 
for  a means  to  minimize  this  loss  function,  the  authors  believe  the 
Graphics  target  clutter  frame-to-f rame  comparison  model  coupled  with 


the  two-dimensional  gradient  and  spatial  discriminator  model  offer  the 
most  effective  and  sophisticated  approach  to  optimize  the  true  target 
selection  criteria  for  most  clutter  conditions. 

However,  it  is  anticipated  that  even  after  maximizing  the  use  of 
these  two  computer  models  in  complex  clutter  environments,  some  loss 
function  may  still  exist.  Consequently,  additional  independent  samples 
of  data  should  be  utilized  with  multispectral  data  reduction  algorithms 
for  further  reduction  of  the  previously  mentioned  loss  function. 
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PRINTV  4 PRTHTX 
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SECURE  2 SHIFT  HO  TTPE 


PROSCAN  P2171  74/7*  0PT»1  TRACE  FTN  4.7*7*355  (9/79/76  17.96. IS.  PACE 


STATISTICS 

PROCHATTIERSTW 


SUBROUTINE  ROHOOE  74/74  OPT«l  TRACE  ETN  4.2.74J55  A5/25/7A  12.5A.il.  PACE 


STATEMENT  labels 


SUBROUTINE  PRTRTX  7h/7*  OPT*l  TRACE  FTN  4.2*74355  09/29/76  12.50.34.  PACE 
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COMMON  BLOCKS  LENGTH 
IPttTX  43 

/ / 10115 


26*1  OHMMA  t6SIVI*<6*S>31iai  ‘ niHHlIVI*  <8)  * N0ISN3MI0 


69 


20  CALL  CHSIZW  (2*21  THERMO  1365 

CALL  RITST#<ie.iS,TAOLi¥l  THERMO  i»*6 

CALL  RITE2tflA50tL20*1023«40«l«ll»*ltl 1HPREPARE C 3T.M)  THERMO  1JL7 

CAtt  RITE2»«35«,Ht#it«*3,0#fl,Jl*-ll»tHrR€€HAM  HATHEHAIICAL  tAOOR*  THERMO  1J66 

lTORf.NI  THERMO  1J69 


S'JflROUTINE  SECURE  7b/7b  OPT«t  TRACE  FT* 
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COHHON  IdEFOR (7) . I AFTER (7) .IT  OP. I PICT (1010  01  THERHO  14JS 


SUdROUf t N€  GETPIC  74/74  OPT*l  TRACE  FTW  4.2*74355  09/29/74 12.S0.42. PAGE 


IOTYP  = SHIET<IMIM)XI ,-141 .AND. 38  THFRHO  1547 

tlN£NT«SMiFTUN<l*{m,-2»T,A**0,J779  THER90  1599 


iaEFOR«.»*IGCOOEUVIO*l»  IHERHO  1687 
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»tHTM«*aS«»RISHT-Vt€f  T> THERHO  17*6 
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JAN  78  G E GOWINS'  H T JACKSON 
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Appendix  C.  LISTING  FOR  PROGRAM  KPLOT 
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»»OS»AN  K PLOT  7L77  L OPT  * 1 PTM  L.  7*71.355  11731777  15. LI. 57.  PAGE 
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I 


I 
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